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Abstract. Simple causal arguments put forward by Kibble and Zurek suggest that the scaling 
behaviour of condensed matter at continuous transitions is related to the familiar universality 
classes of the systems at quasi-equilibrium. Although proposed 25 years ago or more, it is only in 
the last few years that it has been possible to devise experiments from which scaling exponents 
can be determined and in which this scenario can be tested. 

In Ref. [T], an unusually high Kibble-Zurek scaling exponent was reported for spontaneous 
fluxoid production in a single isolated superconducting Nb loop, albeit with low density. Using 
analytic approximations backed up by Langevin simulations, we argue that densities as small 
as these are too low to be attributable to scaling, and are conditioned by the small size of the 
loop. We also reflect on the physical differences between slow quenches and small rings, and 
derive some criteria for these differences, noting that recent work on slow quenches does not 
adequately explain the anomalous behaviour seen here. 



1. Overview 

Causality imposes strong constraints on the non-equilibrium behaviour of condensed matter 
systems, particularly in the vicinity of continuous transitions where, adiabatically, correlation 
lengths become arbitrarily large. Since, in reality, disorder-order transitions take place in a finite 
time, correlation lengths must remain finite; causality imposes a maximum speed at which the 
system can order itself. The greater the transition rate, the more constrained the correlation 
lengths will be. Should the ordered system be frustrated, this frustration will be relieved by the 
spontaneous creation of topological defects, whose separation will reflect this finiteness. 

The suggestion that causality alone could lead to the spontaneous creation of topological 
defects was originally made by Kibble [21 [3], not for condensed matter systems, but for the 
Grand Unified Theories (GUTs) present in the very early universe, which possess a cascade of 
transitions, almost all of which lead to the creation of cosmic strings (field vortices). Zurek, 
independently, had arrived at the same conclusions for the spontaneous formation of defects 
for condensed matter systems [HEIE], and this mechanism is known as the Kibble-Zurek (KZ) 
scenario. This has led to an active collaboration between cosmologists and condensed matter 
physicists over many years, most notably over the classification of defects and the similarities 
in the symmetry-breaking patterns and defects of liquid ^He and GUTs [?]• However, the 
more important aspect of the KZ scenario lies in the prediction of defect densities after a 
transition, since this should be amenable to experimental confirmation in condensed matter 
physics. Unfortunately, we cannot yet say the same for the defects of the early universe. 



The simplest transitions are for systems cooled from one phase to another. Consider a 
temperature quench with quench time tq = Tc{dT / dt)^^ in the vicinity of a continuous thermal 
transition with critical temperature Tc- It is predicted that, for a uniform quench applied to 
a homogeneous and isotropic system, the initial separation ^ between defects takes the scaling 
behaviour 



where and tq are fiducial distance and relaxation time scales, which depend upon the system. 
The important point is that the exponent a is derived from the adiabatic critical exponents and 
therefore shows universality. 

There are several ways to derive ([T|) from causal bounds O [B]. Most simply, we assume 
that the true correlation length tracks the adiabatic correlation length until the latter's rate of 
change is forced to try to exceed the relevant causal speed, whereupon the system freezes in, to 
unfreeze when the correlation length is again changing sufficiently slowly to be causally valid. 
This behaviour has been observed in a variety of simulations (e.g. see Refs. [3 [9]). 

There have been many experimental tests of ([T|) that give results commensurate with it 
without being capable of determining a, although not necessarily for temperature quenches. 
However, some experiments have been able to show scaling behaviour with quench time, in 
particular experiments on Rayleigh-Bernard convection patterns [10], charge density waves in 
TbTes [11] and fluxons in annular Josephson junctions \12\ I13j . Several proposals for future 
experiments have been made, such as experiments for the formation of kink defects in linear 
ionic crystals [H] and quantum transitions in cold gases (e.g. Refs. [T5l [TBI [T7]). 

Although this was the original motive, attempts to see these results as reflections of transitions 
in the early universe have to be taken with extreme caution. The early universe is lacking in 
boundaries and homogeneous and isotropic, to which ([1]) applies. Yet none of these attributes 
are possessed by laboratory condensed matter systems. In particular, inhomogeneity goes hand 
in hand with finite size. Inhomogeneity comes in two kinds. Firstly, that due directly to the 
finiteness of the system e.g. a trapped condensate has variable density |18j . The second is 
concerned with the inhomogeneous application of the quench. It is possible to incorporate such 
behaviours within scaling laws, but the end results are complicated \18\ [19]. 

One way to try to avoid the worst of these problems is to use annular systems, originally 
proposed by Zurek [6]. Such annuli are achievable in superconductors, Josephson junctions and 
in annular traps for cold bosonic gases [T7]. If the circumference of the ring can be taken to be 
larger than its other dimensions for the purposes of domain formation, then the periodicity of 
the ring means we can have a homogeneous system without boundaries. Further, if the ring can 
be kept small, the problem of quench inhomogeneity becomes less pressing. The disadvantage 
of finding less than one defect, on average, is mitigated by ease of observation and the absence 
of defect annihilation. There is yet a further reason for choosing small rings of superconductors. 
The KZ mechanism is not the only reason for the spontaneous generation of magnetic flux in 
superconductors. The long wavelength modes of the magnetic field can freeze in spontaneously, 
almost without regard to the superconductor itself. This mechanism, the Hindmarsh-Rajantie 
[20^ [2T] mechanism gives contributions to the spontaneous flux which have a totally different 
behaviour from ([1]). The effect can be ignored for small superconductors. It was with all these in 
mind that the experiments on fluxon formation in annular Josephson junctions were performed 
[12^ 113] , which showed remarkable allometric behaviour of the form given in ([1]) . 

The anticipated scaling exponent for annular Josephson tunnel junctions was unclear, because 
of ambiguities in the fabrication that may or may not have allowed for proximity effects |22j . 
However, single superconducting annuli should be less ambiguous and a recent experiment on 
low-temperature superconducting annuli showed a higher than expected scaling exponent [IJ. 
This experiment is being revisited to look for possible systematic errors, but it is timely to see 
what effects due to size we might expect in superconducting systems that are smaller than the 
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correlation length at the time of defect formation. We have some understanding of this for 
the quantum transitions of Bose gases \16\ [T7] , but these are controlled by the Gross-Pitaevskii 
equation whereas for superconductors phase ordering is associated with diffusion. 

The purpose of this paper is to advance numerical and analytical arguments that the 
behaviour of an annular system changes from KZ scaling to exponential suppression for defect 
densities as small as those reported in Ref. [Ij. Such densities have not been widely studied in the 
literature, and certainly have not been widely explored for annular systems where the qualitative 
behaviour might be expected to be rather different from that for Bose-Einstein condensation 
experiments. 

In Section [2] we estimate the effect of small size by making a Gaussian approximation for a 
diffusive system on a one-dimensional ring. It suggests an interpolation between the canonical 
scaling behaviour for large rings and exponential damping. This is followed in Section [3] by a 
demonstration that this is confirmed by a Langevin simulation for a simple complex field on 
the ring. More realistic simulations are proposed in Section [H Our conclusions can be found in 
Section [5l 

2. Gaussian approximation 

Although causality does impose bounds, the KZ argument as presented above is too simple. 
It makes more sense to think of defect formation as due to the growth of long-wavelength 
unstable modes, where the growth of these is causally constrained. Surprisingly, treating the 
order parameter field as a Gaussian variable satisfying (short-time) linearised equations, which 
embodies this picture, recreates the scaling laws ^ [23i |24j and some key aspects of the nature 
of the transition [25]. It is borne out not only by the one-dimensional numerical simulations 
presented here, but also by more realistic simulations including the gauge field [26] . 

With due caution, we look to the effect of finite size on the winding number of a complex 
field = e*^^^^ confined to a loop parameterised by the coordinate x. We see this as 
a proxy for a Cooper pair field confined to a small superconducting annulus of circumference 
C. On quenching the system from the normal to the superconducting phase, we can define a 
winding number density n(x) = dxO{x)/2iT, where < x < C is the distance along the annulus. 
The total winding number n around the loop is zero on average, but will have non-zero variance 
(An)^ = (n^), and it is this that we measure. The relevant quantity is the two-point correlation 
function (n (x) n (y)) . We assume that the complex order parameter field ^ = + i(\)2)l\f2 is 
a Gaussian variable, whereby <f)\ and 4>2 are, independently, Gaussian fields. 

The correlation function for the winding number density is now determined by the correlation 
function G{x) for the field components, defined by 

{<^a{x)4>,,{y))=ba^G{\x-y\). (2) 

A simple approximation gives 

(3) 

where / (x) = G(x)/G(0) and (0„ (x) 4 (y)) = ^.^G (|x - y\). 

For strongly coupled dissipative systems the linearised time-dependent Landau-Ginzburg 
equations suggest that, for long rings {C/^ ^ 1) / (r) can be approximated well by the Gaussian 
form / (r) ~ exp(— r^/2^^) where S^r ~ C- Although this expression does not show periodicity, 
for long rings this should not matter. Substituting in ([3]) gives (n^) = 0{C /^r) in the case where 
C is very large and the periodicity in the field can be ignored. This linear dependence of (re^) 
on C is just as we would expect from a random walk in phase along the ring and is assumed as 
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input in the KZ picture. Linearising the equations of motion for fast quenches then gives ([T]) 
with a = 1/4, in agreement with the KZ prediction. 

On the other hand, for short rings (C/^ <C 1) we must have f{x) periodic in x (mod C). To 
implement this we need to replace f{x) by its Jacobi 'd function generalisation 

whence {in?) = C'(e~^'^^^^/'-'^). Thus, rather than the power falloff for large loops we have 
exponential damping 

ln(n2) ^ -A'n'^il/C'^ + const. (5) 

This is in contrast to the proposal made in that the slope be doubled, which is more likely 
for a disc. We should not push this approximation too far but it suggests exponential damping 
in the quench time whenever (n^) <^ 1, whether this be for very slow quenches in large rings or 
faster quenches in small rings. 



3. Numerical simulations 

Our Langevin simulations are for the U(l) scalar field theory described above in 1 + 1 dimensions 
with periodic boundary conditions, on a ring of circumference C. The potential is 

vm') = \m^? + \m\ (6) 

where a{t) is decreased linearly in time as a{t) = —t/Tq {t < tq) and a{t) = — 1 [t > tq) 
to model a slow quench. We start at t = —2tq and continue until t = 4tq, by which time the 
defects have frozen out. 

To study this time-dependent system with dissipation and coupled to a heat bath, we consider 
the second-order Langevin equation 

dV 

dfcPa -dl<Pa + Vdt4>a + ^ = Ca , (7) 

with Gaussian noise 



{Ca{x',t'Kb{x,t)) = 2rjT6{x' -x)6{t' -t)6ab, {Ca(x,t)) = 0. (8) 

We take rj = 1 (so that the system is heavily damped) and T = 0.01. It can be shown that 
such a system satisfies a fluctuation-dissipation theorem. We use a stochastic leapfrog method 
to evolve the equations of motion consistently [27j . 

In this model, there is no dynamical gauge field, a reasonable approximation for the present 
system; we are interested in small rings or slow quenches when relatively few defects form so 
(v?) <C 1. Recent work [28] suggests that the KZ scaling result ([1]) should be modified for long 
times after slow quenches where defect annihilation is of importance. However, as we are dealing 
with flux trapping at very low densities, annihilation dynamics should be unimportant. 

Applying periodic boundary conditions, we the phase angle is obtained through 0{x) = 
Prarg(02(2^) + with a discretised winding number density n{x) = {9{x + 5x) — 

6{x))/Sx, mod 2tt. The main bias in statistics of the winding number obtained from simulations 
will be a slight undercounting due to discretisation effects. One can account for this by checking 
that the results are robust to changes in lattice spacing. 
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Figure 1. Results of simulations, showing change between KZ and exponential suppression 
regimes. The doubly logarithmic plot at left shows KZ scaling with a ~ 0.25 giving way to an 
exponential suppression regime below (n^) ^ 0.2. The same data are plotted again on log-linear 
axes in the plot at right. 



Our numerical results are ensemble averages of 2000 separate simulations for each combination 
of parameters. Typically, the lattice spacing is 6x = 0.125 and the timestep 6t = 0.0125. The 
numbers of lattice sites used is in the region of 192 — 512, giving circumferences C between 24 and 
256 in appropriate units. Qualitative checks of the system's sensitivity to discretisation effects 
were carried out, and decreasing the lattice spacing or the timestep by a factor of 2 did not affect 
the results presented here. The Gaussian noise was generated using the 'ranlux' generator with 
Box-Miiller transformations. In simulating the system, we need to quench the system rather 
slowly (potentially leading to growth of errors from the discretisation of the timestep), while 
keeping it fairly small (motivating a reduction in the lattice spacing). We think that our choice 
of parameters offers a reasonable trade-off. 

In Figure [1] we see behaviour similar to that believed to have been seen in experiment [T]: 
the KZ scaling with a ~ 0.25 gives way at small defect densities, (n?) ~ 0.2, to exponential 
suppression. This is consistent with ([5]). It is easy to see how, for a limited set of data in 
the region of (n^) ~ 0.2, finite- volume effects might be mistaken for a doubling of the scaling 
exponent. 

The behaviour we see in Figure[T]does not particularly support the modified scaling behaviour 
at late times that was recently proposed [28]. In any case, it be unlikely that we would see defect 
interactions substantially alter a texture- forming quench of this kind at such low values of (n^). 

While in experimental settings it is only possible to measure (n?), for the numerical simulation 
we can measure the probability of trapping exactly q fluxons or antifluxons fg. At very small 
trapping probabilities where /g ~ for q>2, then {n?) ~ fi. Figure [2] demonstrates the effect 
of system size on fg, and also on (n^). The results demonstrate that keeping the quench time 
relatively short and simply varying the circumference also supports the Gaussian approximation. 
In other words, slow quenches and small systems show the same exponential suppression of flux 
trapping, as follows from 

4. More realistic simulations: the system in its environment 

More realistic simulations require the incorporation of a dynamical gauge field. It is not clear 
how valid the arguments advanced in the previous sections will be for the case of a thick ring, or 
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Figure 2. At left is a plot of the probability fq of finding a given number of fluxons with system 
size with tq = 24; at right this has been combined into a measurement of {n?). The exponential 
scaling with is clear for small volumes. Our simulations suggest that small volumes and 
slow quenches have the same effect. 



where the system is still annular but is no longer circular. Genuine experiments also typically 
have stray magnetic fields that cannot be truly eliminated. None of these effects can be studied 
in 1 + 1-dimensional simulation as there is no dynamical magnetic field. Simulations of the 
noncompact Abelian Higgs model in 2 + 1 dimensions may be little better [2D]. For these 
reasons we favour treating the system as a thin ring in a larger simulational box containing only 
the U(l) field of vacuum electrodynamics [26]. 

There are no great technical obstacles to such a simulation. Performing a Langevin 
simulation while maintaining the Gauss constraint (having imposed temporal gauge Aq = 0) 
is straightforward and has been used in previous Langevin simulations of similar systems, albeit 
with greater symmetry [29] . 

Indeed, the proposed simulation is not particularly elegant. In checking how external fields 
bias the trapping of flux, the penetration depth introduces a length scale which would be ignored 
for the simple case of a thin film. It is neatest to use periodic boundary conditions, in which we 
leave a gap between the edge of the annulus and the sides of the simulational box; the net flux 
through the system must vanish and so trapped flux must return around the edges of the ring. 

We can measure flux trapping in such a simulation by calculating the magnetic flux through 
the ring or by measuring the scalar fleld winding around the inner boundary of the system. The 
results from both methods are in agreement based on preliminary results [26j. 

Perhaps, if we wish to examine finite-size issues in effective models of planar (non- annular) 
superconducting films, the most elegant way may be to create a film in the same way [30]. Then 
a direct parallel can be drawn with the experiments carried out in Bose-Einstein condensates 
[161 [T7] . Periodic boundary conditions may be employed without hesitation so long as a gap is 
left at the edges of the fllm. 

5. Conclusions 

We have performed numerical simulations for small annular systems with broken U{1) symmetry, 
an idealisation of superconducting rings. We have demonstrated that the simple power behaviour 
of ([1]) breaks down when the probability of flnding a single defect falls below about 0.2, more 
or less independently of the microscopic properties of the model. As we anticipated in a simple 



Gaussian approximation the likelihood of finding defects is then seen to become exponentially 
suppressed with quench rate or system size. This behaviour is very different from that of defect 
production in small Bosc-Einstein condensates. 

We have also discussed a potentially more valid simulational protocol for more realistic 
(planar) annular superconductors, incorporating the three dimensional electromagnetic field 
(and an external applied field) in a consistent way, and explained some of the limitations of such 
an approach. 
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